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i 

A  population  of  susceptible  individuals  exposed  to  m  contagious  diseases 
is  considered.  The  progress  of  this  epidemic  among  the  individuals  is  modeled 
by  an  m-dimensional  stochastic  process.  The  components  of  this  process  repre¬ 
sent  the  number  of  infective  individuals  with  the  respective  diseases  at  time  t . 

A  class  of  m-dimensional  stochastic  processes  is  constructed.  These 
processes  describe  the  progress  of  the  epidemic  models  considered  in  the  sequel . 
Exact  and  approximate  formulas  for  the  joint  and  marginal  state  probabilities 
of  these  models  are  obtained.  It  is  shewn  that  the  approximate  formulas  are 
relatively  simple  functions  of  time  while,  the  derivations  of  the  exact  for¬ 
mulas  involve  tedious  computations.  The  results  obtained  in  the  paper  are 
applied  to  a  sample  of  examples.  v 

Key  words:  m-dimensional  simple  epidemics,  exponential,  negative  bino¬ 
mial  and  Poisson  random  variables,  stochastic  processes  and  convergence  in 
distribution. 
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1 .  Introduction  And  Summary . 

In  a  simple  epidemic  situation  we  assume  that  a  population  of  suscepti¬ 
ble  individuals  (susceptibles)  is  exposed  only  to  one  contagious  disease 
(disease)  [Bailey  (1975)].  However,  frequently  susceptibles  are  exposed 
simultaneously  to  more  than  one  disease,  as  is  the  case  with  different 
types  of  flu.  In  this  paper  we  consider  a  population  of  susceptibles 
exposed  to  m  diseases,  m  =  1,  2,  ...  .  We  say  that  this  population  undergoes 
an  m-dimensional  simple  epidemic  if  the  following  six  assumptions  hold. 

(1.1)  At  each  point  in  time  at  most  one  susceptible  contracts  a  disease. 

(1.2)  Each  susceptible  constracts  at  most  one  disease. 

(1.3)  Once  a  susceptible  contracts  disease  r,  r  =  1,  ....  in,  he  remains 
contagious  for  the  duration  of  the  epidemic. 

(1.4)  An  infective  individual  (infective)  with  disease  r,  r  *  1,  m,  can 

only  transmit  that  disease. 

(1.5)  All  interactions  between  a  susceptible  and  an  infective  with  a  specified 
disease  are  equally  likely  to  result  in  an  infection. 

(1.6)  Individuals  neither  join  nor  do  they  depart  from  the  population. 

Note  that  for  m  =  1  the  m-dimensional  epidemic  reduces  to  the  univariate 
simple  epidemic. 

Let  Tq  denote  the  first  time  we  have  at  least  one  infective  with  each 
of  the  various  diseases,  and  n,  finite  or  infinite,  the  number  of  suscepti¬ 
bles  at  time  Tq.  We  describe  the  progress  of  m  dimensional  simple  epidemics 
among  susceptibles  by  an  m-dimemsional  stochastic  process 
X^t)  ■  CXn  j(t),...,Xn  m(t)],  t«[0,»).  The  components  of  X^(t)  represent 
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the  number  of  infectives  with  the  respective  diseases  at  time  t  measured 
from  Tg.  In  Section  2  we  construct  a  variety  of  m-dimensional  stochastic 
processes.  These  processes  describe  the  progress  of  the  epidemic  models 
considered  in  the  sequel. 

The  computation  of  state  probabilities  in  epidemic  models  is  of  major 

interest  to  researchers.  In  Section  3  we  derive  formulas  for  the  joint  and 

marginal  state  probabilities:  P  ,  (t)  =  P{X  (t)-X  (0)=k>,  and 

n,K  “H 


Pn,k,r(t)  =  PtXn,r^t^'Xn,r^°^'=k}*  where  *  =  ,  • . .  ,km3,  k,  kj,  ...»  kme{0,l, 

r  ■  1,  . ...  m,  and  te(0,“).  We  calculate  these  formulas  without  the  tradi¬ 
tional  use  of  the  differential  equations  associated  with  the  state  probabil¬ 
ities.  This  is  done  by  utilizing  a  formula  for  the  distribution  function 
(df)  of  a  sum  of  independent  exponential  random  variables  (rva's)  given  by 
Billard,  Lacayo  and  Langberg  (BLL)  (1980). 

The  for  :i.  .;  for  the  state  probabilities  obtained  in  Section  3  are 
rather  complicated.  To  overcome  this  deficiency  we  derive  relatively  simple 
approximations  to  the  state  probabilities  when  the  initial  number  of  suscep- 
tibles:  n,  is  sufficiently  large.  To  be  more  specific  some  notation  and  a 
definition  are  needed.  Let  X  (t),  X  (t),  be  m-dimensional  stochastic 

processes,  X^Ct)  -  Xn(t)  "  ^r-1  Xn,r(t)’  the  total  number  of 

infectives  at  time  t  measured  from  TQ,  and  let  Y(t)  be  a  nonnegative  rv, 
tcC0,»),  n  ■  1,  2,  ...  . 

Definition  1.1.  Let  n  ■  1,  2,  Then  the  transition  rate  of 

disease  r,  r  «  1,  ...»  m,  at  time  t,  tc[0,*»)  is  given  by 

*im  ft*h)  -  X„  ft)  •  l|X  (t)> ,  and  is  denoted  by  R(X  (t),r). 

“>*  **|*  **  *■ 
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In  Section  4  we  assume  that  for  all  k^  =  [k^,...k  ],  an<*  given  s,  6e  (0,®) 

(1.7)  am  X  (0)  =  XJO), 

n->® 

(1.8)  am  RQC^O^k.r)  =  R(X^(0)+k  ,r) ,  r  *  1,  ....  m, 
n-*» 

(1.9)  E(Y(s)}6  <  ®,  and  that 


(1.10)  for  xeCX^CO) ,«)  and  almost  all  ne(l,2,...},  P{Xn(s)2x)  £  P(Y(s)ax}. 

We  show  that  under  Conditions  (1.7),  (1.8),  X^(t)  converges  in  distribution  as 

n  +  •  to  X  (t)  for  all  te(0,®).  Using  this  result  we  can  approximate  P  ,  (t), 

and  P  .  (t)  by  the  respective  state  probabilities  of  the  process  X  (t).  Further, 

n,K,r 

we  show  that  £im  E{X  (s)}6  =  E(X  (s)}S,  and  that  Urn  E(X  (s))P  =  E(X  (s)}6, 

n  °°  _  n,r  00  ,r 

n-H*>  rt+®  ’ 

r  *  1,  ....  m,  provided  Conditions  (1.7)  through  (1.10)  hold.  Note  that  if 
6eC2,“),  we  can  approximate  EX^  r(s),  EX^(s),  VariX^  r(s)},  and  Var{Xn(s)},  by 
the  equivalent  counterparts  of  X^s). 

BLL(1979)  consider  a  special  class  of  m-dimensional  simple  epidemic  models 
and  name  them  the  symmetric  m-dimensional  simple  epidemics.  We  say  that  a  pop¬ 
ulation  exposed  to  m  diseases  undergoes  a  symmetric  m-dimensional  simple  epidemic 
if  Assumptions  (1.1)  through  (1.6)  hold  and  if  the  transition  rate  of  disease  r, 
r  »  1,  ...,  m,  at  time  t,  t€[0,<*>)  is  given  by: 


(1.11)  R(X„(t),r)  = 


n"laXn,r(t)(n+Xn,r(0)'Xn,r(t))’  n  * 


aX  (t) 
®  ,r 


n  = 


where  o€(0,«).  The  transition  rates  given  in  Equation  (1.11)  for  the  case  n  =  ® 
depend  only  on  the  number  of  infectives  with  the  respective  diseases.  This 
special  structure  of  the  transition  rates  is  shared  by  a  variety  of  models  as 
is  shown  in  the  last  section.  In  Section  5  we  assume  that 
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(1.12)  R(X.(t),r)  =  A (X^  (t) ,r) ,  r  =  1,  ....  m,  te[0,»), 

r- 

where  A(k,r),  k  =  x _ (0),  X_  _(0)  +  1,  ....  r  =  1 . m,  are  sequences 

>r  |P 

of  positive  real  numbers.  We  prove  that 

(1.13)  for  r  =  1,  ....  m,  k  =  0,  1,  ....  and  te(0,“>), 

P{X-fr(t)-X.tP(‘>)^k}  = 

=p^5=lCA(X.,r(0)+c»-1'r5j"1Vt} 

where  U^.U^,  ....  are  independent  indentically  distributed  (iid)  exponential 
rva's  with  mean  1.  Further  we  show,  as  expected,  that 

(1.14)  for  all  te(0,“),  the  components  of  ^(t)  are  independent. 

These  results  are  applied  in  Section  6  to  simplify  the  approximations  to  the 
state  probabilities  and  to  related  moments.  Note  that  under  Condition  (1.12) 
the  processes  X^  ^(t),  . . . ,  X^  m(t),  describe  the  progress  of  m  independent 
univariate  simple  epidemic  models  among  susceptibles  in  closed  populations 
of  infinite  sizes. 

In  Section  6  we  apply  the  theory  developed  in  the  previous  sections  to 
some  classes  of  m-dimensional  simple  epidemics,  and  other  models  in  Epidemiology. 


*ir- 


1 
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2.  Model  Construction. 

In  this  section  we  construct  m-dimensional  stochastic  processes.  These 
processes  describe  the  progress  of  m-dimensional  simple  epidemics  in  a  pop¬ 
ulation  of  a  finite  or  infinite  number  of  susceptibles . 

First,  we  introduce  some  notation.  Let  N(n)  ■  n  +  1  for  n  <  «,  and 
-  •  for  n  *  «,  and  T.  1  £  k  <  N(n) ,  the  k—  interinfection  time  defined  as 
the  time  that  elapses  between  the  X  (0)  ♦  k  -  1  and  the  Xn(0)  ♦  k  infection. 
Further,  let  £  ,  ,  1  £  k  <  N(n) ,  be  a  rva  assuming  values  in  the  set  {l,...,m}, 

and  £  _  =  0.  The  rva  £  .,  1  £  k  <  N(n),  designates  the  disease  responsible 
ti  i  u  n  ^ic 

for  the  Xn(0)  *  k  infection.  Finally,  let  I_ be  the  indicator  function, 
rk 

^n,0  *  Sn,k  "  ^q=l  Tn,q*  and  Sn,N(n)  *  ».  The  rva  S  .  ,  1  £  k  <  N(n), 

n  i  a 

describes  the  time  measured  from  TQ  until  the  *n(0)  +  k  infection.  In  par- 
represents  the  duration  time  of  the  m-dimensional  simple 

epidemic. 

Note  that  for  0  £  k  <  N(n) ,  r  *  1 ,  . ..,  m,  and  te(0,«),  the  following 
event  equality  holds. 


(2.1) 


(X  (t)-X  (0)=k) 

v  n,rv  *  n,r v  ’ 


U  (Sn  qSt<Sn  a*l'  Ej-lI(5n  j"r>k> 
k£q<N (n)  n,q  n,q  1  3  1  n,} 


Thus,  to  construct  the  process  X^(t)  if  suffices  to  determine  the  df 's  of  the 

random  vectors  (rve's)  [T_  _,£_  _,q«l, . . . ,k),  1  £  k  <  N(n).  Next  we  determi... 
-■■■  ■  *  n,q  n,q  •  * 

the  df's  of  these  rve's.  We  need  some  more  notation.  Let  J  .  be  the 

jnjLtl 

index  set  of  all  infectives  with  disease  r  at  time  TQ  ♦  SR  and 

c.  ■  Xn  r(0)  ♦  ^”gI(Cn,q“r)  *  1  £  k  <  N(n) ,  r  ■  1,  ...»  a.  The  rva  Cn^r 
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represents  the  number  of  infectives  with  disease  r  at  time  +  S  ,  , .  Fur- 

0  n,k-l 

ther,  let  u(n,k,q,r) .  1  s  k,q  <  N(n),  r  =  1,  . ..,  m,  be  positive  real  numbers 

and  t  ,  .,  i  =>  1,  ....  X  (0)  +  k  -  1 ,  1  s  k  <  N(n),  be  rva's  defined  on  a 
n  p  x  1 1\  u 

common  probability  space.  The  rva  t  .  .  describes  the  time  measured  from 

n,i,K 

the  X^CO)  +  k  -  1  infection  until  the  i—  contagious  individual  causes  the 
Xn(0)  +  k  infection.  Throughout  we  assume  that 


(2.2) 


the  conditional  rva',  (rn>1>k | . 

i  =  1,  ...i  X  (0)  ♦  k  -  1,  1  s  k  <  N(n),  are 
independent  exponentially  distributed,  and  that 


(2.3)  E{Tn>i,k^n,0***"5n,k-l}  = 


-1 


"  u  (n.k,Cn>k>r,r),  i^n>k<r»  1^<  N(n) 


1,  • . « ,  m . 


We  are  ready  now  to  determine  the  df's  of  the  rve's  [Tn  ^,En  ^,q  =l,...,k), 

1  s  k  <  N(n).  First,  we  determine  the  df's  of  the  rve's  [£  . ], 

n  f  x  n  f  k 

1  s  k  <  N(n),  We  need  the  following. 

Lemma  2.1.  Let  us  assume  that  Conditions  (2.2),  (2.3)  hold,  let 
r  *  1,  ....  m,  and  1  s  k  <  N(n).  Then 

(2.4)  ■,U„,k-*l«n,0.-.«n>k.1>  - 

’*'n,k,rll^n'lc,Cn,kjT'r^^£«lCn,k,£u^n’k'Cn,k,£’^  ^  ' 

Proof.  Note  that  the  event  (En  k*r)  is  equal  to 

*  m 

(mintt  .  .  :  i«J  .  ,}<min{T  .  .  :  ie  U  J  ,  „}).  Consequently  Statement  (2.4) 
n,i,K  n,x,r  n,i,K 


follows  by  Conditions  (2.2),  (2.3)  and  simple  integral  evaluations.  || 

Note  that  Equation  (2.4)  determines  the  df's  of  the  rve's  [E  E_  .], 

n | i  n,K 

1  l  k  <  N(n).  Next,  we  determine  the  df's  of  the  conditional  rve's 
^Tn,q^n,o*"*^n,q-l*'  q  "  1  s  k  <  N(n). 


Lemnr  2.2.  Let  us  assume  that  Conditions  (2.2),  (2.3)  hold.  Then 

the  conditional  rva's  (T  ,  |£  n,...,S  ,  } ,  1  <  k  <  N(n),  are  indepen- 

n,K  n,u  n  j  k  ~ l 

dent  exponentially  distributed  with  means  equal  respectively  to 
C^r=lCn,k,rp (n'k,Cn,k ,r'r) ]  *  * 

Proof.  Note  that  T  ,  -  mir.fr  .  ,  :  i  =  1 . X  (0)+k-l>,  Is  k  <  N(n) 

-  n,k  n,i,k  '  n 

Consequently  the  result  of  the  luma  follows  by  Conditions  (2.2),  (2.3)  and 
some  simple  properties  of  exponential  rva's.  || 

Clearly  Lemmas  2.1  ar.d  2.2  together  determine  the  df's  of  the  rve's 
[T  o,Sn  1  S  k  <  N(n). 

Finally,  wc  show  that  the  transition  rates  of  an  m-dimensional  simple 
epidemic,  as  ejected,  determine  uniquely  the  epidemic.  We  need  the  follow! 
lemma. 

Lemma  2.5.  Let  ur  assume  that  Conditions  (2.2),  (2.3)  hold.  Then 
RQ^OO.r)  =  Xnjr(t)p(n,Xn(t),Xn  r(t),r),  r  =  1,  ....  m,  te[0,cc). 

Proof.  The  result  of  the  lemma  f  .  by  the  mcmoryless  property  of 
exponential  rva's  [Bariev/,  Pvoschan  (1975),  p.  56];  Equation  (2.1)  and 
Conditions  (2.2),  (2.3).  || 

Consequently  oy  Lemr.:a  2.3  and  Condition  (2.3)  we  obtain  for  ieJ  , 

n,K ,r , 

1  i  k  <  H(n) ,  r  =  1,  ...,  m,  that 


(2.5) 


E{Tn,i,k^n,0''”’Cn,k-l}  = 

3  Cr.,k,rCR({Cr.,k,l’--*»Cn,k,m}  J 


Thus,  clearly  m-dimcnsional  simple  epidemic  models  are  determined  by  their 
transition  rates. 


f 


s 


3.  Formulas  for  the  State  Probabilities. 


Let  k,  kj,  ....  kme{0,...,n},  n  =  1,  2,  ....  k  -  . ,k^] , 

s(k)  =  £™_jk  ,  and  6^,  6^,  ....  be  positive  real  numbers. 

This  section  contains  formulas  for  the  joint  and  marginal  state  prob¬ 
abilities:  P  ,  (t) ,  and  P  .  (t),  r  =  1,  . . . ,  m,  te (0 ,<“) .  These  formulas 

n,K  n jK  ,r 

are  calculated  without  the  traditional  use  of  the  differential  equations 

associated  with  the  state  probabilities.  Rather,  we  utilize  an  available 

formula  for  the  distribution  function  of  a  sum  of  independent  exponential 

rva's.  For  the  sake  of  completeness  we  present  this  formula. 

Theorem  3.1.  [BLL  (1980),  Theorem  23.  Let  M  be  a  positive  integer, 
m  jn+1 

and  -f^(j)  *  J  lie"  ,  j  =  M,  M+l ,  ...  .  Then  for  te(0,°°) 

j1+...+jM=3  <1=1 

To  aid  in  computing  the  joint  and  marginal  state  probabilities  we 
introduce  the  following  notation.  Let  £(k)  =  [£Q, . . .  ,£^3,^=0 ,1^  , . . .  , . .  ,m}  , 

B(k)  *  U(s(k)):^I(£q=r)=kr,r-l . m>,  and  A(r,j,k)  =  tf(j )  :ZjmQI  ^e-r)-k> , 

j  a  k,  k+1,  ...  . 

By  Equation  (2.1)  we  obtain  that  for  te(0,«) 


(3,2)  Pn,k(t)  =  £(.s(|))CP{Sn,s(k)St<Sn,s(k)+llCn,q='eq'q=0',,-,s(y}' 

•P(Cn  q=£j,q=0, . . . ,s (k)}I (£(s (k))eB(k)) ],  and  that  for  r  =  1,  ...,  m 

<3*3>  Pn.k,r^  *  5-k  ^)CP{Sn.j<t<Sn.j*ll«n.e^e--0 . J>‘ 

•PUn>e=*e.e»0,...,j}I(£(j)€A(r,;j,k))3. 

Thus,  to  compute  P  .  (t)  and  P  ,  (t)  it  suffices  to  evaluate  P{£  <*l  ,q=0,..,,k}, 

n,K,r 

.nd  l’(Snik«<Snjk,1|{n(q.iq.,.0 . k>. 
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Now,  we  present  formulas  for  these  probabilities.  Let 

D  .  (£(k),q)  =  X  (0)  +  yq'J,  I(£.=r),  and 

n,k,r  —  2__  n»r  J=0  J 

n(n,k,£(k) ,q)  =  Ir=lu^n,q,nn,k  r(^k)>q)>r)*  q  =  1  *  •••>  k  *  1  •  First,  by 
Statement  (2.4) 


(3.4)  PUn,q=Vq=0’"',k}  = 

k’ 


n  [(n(n,k,£(k) ,q) }  1  n  (u(n,q,D  (£(k),q,r) 
=1  _  r=l  n>K,r 


Hl*r) 


Next,  let  fk(ei,...,ek,t)  be  the  density  function  of  the  rva  Then 

for  k  =  1,  2,  . ...  t  e(0,«) 

0-5)  ' 


=  /  e 


fk i » • • • »9k»u)du  -  0k+1f(ei,...,0k+1,t) 


Thus ,  by  Lemma  2 . 2  and  Theorem  3 . 1 


P{Sn,kSt<Sn,k+l I Cn,q=£q=0’ • • ‘ 'k}  = 
-n(n,k,£(k) ,l)t 


k  =  0 


n.n^ — — L  +  +j  n  {n(n,k,£(k),q)}jq+1  k  *  n 

J~n  j*  31  q=l 

(n(n>k,f(k),ktl)}~1rk  l  +  +.  • 

3  (j-l)!  31  **•  3k+l  J 


k+1  j  +1 

n  {n(n,k,£(k)  ,q)}  q 
q=l 


1  <  k  <  n. 


Consequently,  the  formulas  for  the  joint  and  marginal  state  probabilities 
can  be  obtained  from  (3.2)  and  (3.3) by  substitution.  Note  that,  as  expected, 
the  exact  formulas  for  the  joint  and  marginal  state  probabilities  are  rather 


complicated. 
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4.  The  Asymptotic  Approximations . 


In  this  section  we  present  the  asymptotic  approximations  of  the  joint 
and  marginal  state  probabilities  and  of  some  related  moments.  All  limits  are 
calculated  as  n  -*■  <*>.  First,  we  show  that  under  Conditions  (1.7),  (1.8),  X^(t) 
converges  in  distribution  to  X^t)  for  all  t  e(0,“).  Then  we  prove  that 
£im  E(Xn(s)}e  =  E{Xoo(s)},  and  that  £im  E{Xn  r(s)}6  =■  EiX^  r(s)}e  for  r  *  1,  .... 
provided  Conditions  (1.7)  through  (1.10)  hold. 

To  establish  the  convergence  of  X^(t)  to  X^t)  we  need  the  following  two 
lemmas . 

Lemma  4,1.  Let  us  assume  that  Conditions  (1.7),  (1.8)  hold,  let 
k  »  1,  2,  ...,  and  Z^,  . . . ,  Z^  e{ 1 , . . . ,m) .  Then 

(4.1)  Zim  PU  =£  ,q-l,...,k)  =  P{£„  =Z  q=l,...,k}. 

11,4  4  >4  4 

Proof.  Let  Zn  =  0.  Then  P{£  -Z  ,q=l,...,k}  = 

-  0  n,q  q 


m. 


n  P{E  =Z  Ic  .-Z. ,  j=0, . . . ,q-l)  .  Consequently  Statement  (4.1)  follows  by 

qBl  n»q  q  J 

Conditions  (1.7),  (1.8),  and  Statements  (2.4),  (2.5).  || 

Lemma  4.2.  Let  us  assume  that  Conditions  (1.7),  (1.8)  hold,  and  let 

k  *  1,  2,  ...  .  Then  the  conditional  rva's  (S  ,  |£  ,q=0, . . . ,k-l)  converge 

n  |  k  n  1 4 

in  distribution  to  the  conditional  rva  (S  ,  |c  q=0, . . .  ,k-l)  . 

”>*  °°,q 

Proof.  To  prove  the  result  of  the  lemma  it  suffices,  by  the  Cramer-Wold 

device  [Billingsley  (1968),  p.  49],  to  show  that  the  conditional  rve's 

{ CT  .,..., T  .  ]Je  ,q*0, . . . ,k-l)  converge  in  distribution  to  the  conditional  rve 
n ,  i  n,K  n,q 

(CT^  2,...,Tm  q=^-q»qa0» •  •  •  *k-l) .  The  preceeding  statement  follows  by 

Conditions  (1.7),  (1.8),  Lemma  2.2,  and  Statement  (2.5).  || 
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We  are  ready  now  to  show  that  X^(t)  converges  to  X^tt). 

Theorem  4.3.  Let  us  assume  that  Conditions  (1.7),  (1.8)  hold  and  let 
t  c  (0 , °° ) .  Then  the  rve's  X^tt)  converge  in  distribution  to  the  rve  ^(t). 

Proof.  Note  that  X_(t),  X_(t),  a  =  1,  2,  ....  are  discrete  rve's.  Thus. 

—co  —n 

[Billingsley  (1968),  p.  16],  to  prove  the  result  of  the  theorem  it  suffices  to 
show  that  for  all  1<  =  [k^,...^],  kj,  km  e{0,l,...} 

(4.2)  lim  P{X  (t)-X  (0)=k)  =  P{X  (t)-X  (0)=k). 

Statement  (4.2)  follows  by  Equation  (3.2),  Conditions  (1.7),  (1.8)  and  Lemmas 
4.1,  4.2.  | | 

By  Cramer-Wold  device  and  Theorem  4.3  we  obtain  that 

Corollary  4.4.  Let  us  assume  that  Conditions  (1.7),  (1.8)  hold,  let 
t  e  (0,“>) ,  and  r  =  1,  . ..,  m.  Then  (i)  X  (t)  converges  in  distribution  to 

Ti  ,T 

X^  r(t),  and  (ii)  XR(t)  converges  in  distribution  to  X^t). 

Clearly  Theorem  4.3  and  Corollary  4.4(i)  insure  that  the  joint  and  marginal 
state  probabilities:  Pnk(t),  Pn>k>r(t),  can  be  approximate  by  POT>k(t),  P„  k  r 
respectively. 

Now,  we  establish  the  moments  convergence.  First,  we  show  that 

iim  E{X  (s)}B  =  E{X  (s)}B. 
n  ® 

Theorem  4.5.  Let  us  assume  that  Conditions  (1.7)  through  (1.10)  hold. 

Then  £im  E{X  (s)}B  •  E(X  (s)}B. 
n  ® 

Proof.  Note  that  E{X  (s)}B*B/*yB  *P{X  (s)>y}dy.  Thus,  the  result  of  the 

n  0  n 

theorem  follows  by  Corollary  4.4(ii),  Conditions  (1.9),  (1.10)  and  the 
dominated  convergence  theorem  [Loeve  (1963),  p.  125],  j| 

o 

Next,  we  show  the  convergence  of  E{X  (s)}  . 

n,r 


(t). 
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Theorem  4.6.  Let  us  assume  that  Conditions  (1.7)  through  (1.10)  hold, 

and  let  r  *  1,  ....  m.  Then  £im(X  (s ) } ®  =  E{X  (s)}^. 

n.r1-  1  ®,r'‘  1 

Proof.  Note  that  E{X  (s)}^  =  &f  y®  1P{X  (s)>}dy,  and  that  for 

li ,  r  q  n  >  a 

ye(-®,“),  n  =  1,  2,  ...,  P{Xn  r(s)>y)  s  P{Xn(s)>y).  Thus,  the  result  of  the 
theorem  follows  by  Corollary  4.4(i),  Conditions  (1.9),  (1.10),  and  the 
dominated  convergence  theorem.  |  j 

Note  that  if  Condition  (1.9)  holds  then  E{Y(s)}^<«  for  6e(-«,0J.  Now, 

assume  that  Conditions  (1.9),  (1.10)  hold  for  all  se(0,»)  and  some  6e[2,“). 

Then  we  can  approximate  EX  (t),  EX  (t),  Var{X  (t)},  and  Var{X  (t))  by  the 

n  n  jT  n 

equivalent  counterparts  of  X  (t)  for  all  te(0,«),  and  r  =  1,  ...,  m. 


1  7 


5.  Some  Special  m-bimensional  Simple  i.pidemic  .'.omi'ls  . 

Let  X^ft),  te[0,"),  Lie  an  m-dimensional  stochastic  process  describing 

the  progress  of  an  ra-dimcnsional  simple  epidemic  among  infinitely  many 

susceptibles .  We  assume  that  the  transition  rates  of  X^(t)  depend  only  on 

the  number  of  infectives  with  the  respective  diseases.  Under  this  assumption 

we  identify  the  distribution  functions  of  X  ,(t),  ....  X  ft),  and  show,  as 

°°,1  ’  °°,rav  '  ’  ’ 

expected,  that  X^ft)  has  independent  components  for  te(0,“).  These  results 
are  then  used  in  Section  6  to  obtain  relatively  simple  approximations  to  the 
state  probabilities  and  to  some  related  moments  of  certain  classes  of 
m-dimensional  simple  epidemics. 

We  need  some  notation.  Let  Zf(t),  te[0,°»),  r  =  1,  . ..,  m,  be  independent 

stochastic  processes  assuming  values  in  the  sets  {X^  (0) .X^  (0)+l,...), 

®  ,r  “ »r 

respectively  and  determined  by  the  following  two  equations. 

(5.1)  Zr(0)  =  X^fO),  r  =  1.  ....  b. 

(5.2)  Hzr(t)-^  - 

-  P{I1q=;1[A(Zr(0)+q-l,r)]-1UqSt} 

te(0,«),  r  =  1,  ...,  m,  k  =  1,  2,  ...  . 

Further,  let  Z(t)  =  Z„(t),  W  (t)  =  0,  W  (t)  =  Inf{t:  Z(t+][k'lWo)-Z(0)-k} , 

-1  *  0  -  q*0  q 

t€(0,~),  and  D,.  -  )“  _  W  ,  k  =  0,  1,  ...  .  Finally,  let  X.  ,  k  *  1,  2,  . . . ,  be 

A  CJ  Iv 

rva's  assuming  values  in  the  set  {l,...,m>,  determined  by  the  following  set 
equality 

(5.3)  (Xk=r)  =  (Z£(Dk)  -  Z£(Dk_1)  =  0, 

Z^(Dk)  -  Zr  (Dk_i)  =  1,  r  .  ,m) ,  r  =  1,  .  •*,  m. 


The  processes  Z^t),  Z(t),  describe  the  progress  of  m  independent  simple 

epidemics  among  susceptibles  in  closed  populations  of  infinite  sizes.  Assume 
we  observe  the  progress  of  these  m  independent  simple  epidemics.  Then  Z(t) 
represents  the  total  number  of  infectives  at  time  t,  to'^t)  models  the  time 
that  elapses  between  the  Z(0)  +  k  -  1  and  the  Z(0)  +  k  infection,  describes 
the  time  until  the  Z(0)  +  k  infection,  and  \ ^  indicated  the  simple  epidemic 
responsible  for  the  Z(0)  +  k  infection. 

In  this  section  we  assume  that  Condition  (1.12)  holds,  and  prove  that  the 
processes  X^Ct)  and  Z_(t)  =  [Z^  (t)  , . . .  ,Zm(t)  ] ,  te[0,«>),  are  stochastically  equal. 
In  particular  we  obtain  Statements  (1.13),  (1.14). 

Let  t£(0,®),  k  =  1,  2,  ....  and  r  =  1,  ...,  m.  Then  the  following  set  equality 
holds . 


(5.4) 


(Zr(t)-Zr(0)>k)  = 


=  U,  (D  St<D  ,  I  (A  ,*r)£k)  . 

q=kv  q  q+l’^jal  j 


By  the  set  Equalities  (2.1),  (5.4),  the  processes  X^ft),  Z_(t)  are  determined 

by  [T^  ,5^  ,q=l  ,2, . . . ],  and  [W  ,X  ,q=l,2,...],  respectively.  Thus,  to  prove 

that  X„(t),  Z(t)  are  stochastically  equal  it  suffices  to  show  that 

CT  ,C_  ,q*l,2,...],  [W  1  ,q=l,2,...],  are  stochastically  equal. 

>M  H  R 

Next,  we  prove  the  preceeding  statement.  We  need  the  following. 

Lemma  5.1.  Let  us  assume  that  Condition  (1.12)  holds.  Then  CT  , ,3 
1  •  ®,i  “>,1 

and  CW^ , A j ]  are  equal  in  distribution. 

Proof.  Let  te(0,«)>  and  r  •  1,  ...,  m.  Then  PCW^t.A^r)  * 

-  P{t<CA(X->r(0),r)]'1U1<Cd(Xa>^(0),£)3"1lJt,r^l,...,m}=  P{Tw  j-r).  || 


We  are  ready  now  to  prove  the  main  result  of  this  section. 

Theorem  5.2.  Let  us  assume  that  Condition  (1.12)  holds.  Then  the  two 

sequences  [T^  ^q=l,2,...J  and  [IV A  ,q=l,2,...J,  are  stochastically  equal. 

Proof .  We  have  to  show  that  for  all  k  =  1,  2,  ....  the  eve's 

[Tm  ,5^  ,q=l,...k_i,  [W  ,X  q=l . . . ,kj  are  equal  in  distribution.  We  prove  this 

» 4  »4  4  4 

statement  by  an  induction  urguement  on  k. 

First,  we  present  some  notation.  Let  (  =  1,  ...,  m,  Z(t,£)  =  Z(t+Wj), 
Zr(0,fl)  =  Zr(0),  Ur  =  1,  ....  m,  and  2^0,1)  =  Z^(0)  ♦  1.  Further, 
let  Wk(£),  Xk(£),  be  defined  as  W'k ,  X^,  respectively,  k  =  1,  2,  ...,  where 
Z[t,U  replaces  2^(t)  .  Finally,  let  X^(t  ,£),  te(0,®),  be  a  m-diraensional 
stochastic  process  as  defined  in  Section  2  such  that  X^  (U,t)  =  (0), 

Ur  =  1,  ....  m,  X^^O.-C)  =  X^^CO)  *  1,  and  let  T„>k(f),  k  =  1,  2,  ... 

be  the  interinfection  times  and  infection  causes  associated  with  X^Ct.f). 

We  now  return  to  the  proof  of  the  theorem.  Let  k  =  2,  3,  ...,  and  te(0  ,«■>). 
By  the  memoryless  property  of  exponential  rva's  the  two  conditional  rve's 
{[Ta.>q>5n,jCl»cl=2»*-*»k]|T«>1=t=  and  {[Wq,Xq,q=2,...,k}|wi=t,X1=f},  are 

equal  in  distribution  to  the  two  rve’s  [Tw  (£ ) ,q:l , . . . ,k-l ]  and 

(f)  ,Xq (£)  ,q=l , . . .  ,k-l  ]  respectively.  Consequently  the  result  of  the  theorem 
follows  by  Lemma  5.1,  Bayes  formula,  and  an  induction  arguement  on  k.  || 

For  reference  purposes  we  summarize  the  main  result  of  this  section  in  the 
following  theorem. 

Theorem  5.3.  Let  us  assume  that  the  transition  rates  of  the  m-dimensional 


process  X^Ct),  teCO,®),  are  given  by  Equation  (1.12).  Then  Statements  (1.13), 
(1.14)  hold. 


It 


6 .  Examples . 

Now  we  apply  the  theory  developed  i:i  the  previous  sections  to  some 
specific  m-dimensional  simple  epidemic  classes.  Throughout  we  assume  that 
Condition  (1.7)  holds. 

Following  the  classical  approach  to  the  univariate  simple  epidemic 
[Bailey  (1975)]  one  can  assume  that  the  transition  rates  of  X^(t)  are  propor¬ 
tional  to  the  number  of  infcctives  with  the  respective  diseases  and  the  total 
number  of  susceptibles .  Let  ,  . ..,  a^£(0,“>).  Then  the  transition  rates  of 
the  classical  m-dimensional  simple  epidemic  are  given  by 


(6.1)  RQCJtJ.r)  =  ar  X  (t)  (n-Xn(t)*Xn(0)) ,  r  «  1,  ...,  m  ,  te[0,«). 

From  Lemmas  2.2,  2.3  the  duration  time  of  a  classical  m-dimensional  simple 
epidemic  is  less  than  or  equal  to  [minta^ :r=l , . . . ,m) J  (Xn(0)+q-l) (n-q+1 ) ] 

Thus,  the  duration  time  tends  to  zero  as  n  -*•  ».  Consequently  the  classical 
models  are  not  describing  properly  the  progress  of  the  epidemic  when  n  is 
large.  BLL(1979)  adjust  the  classical  transition  rates  and  use  the  symmetric 
m-dimensional  simple  epidemic  to  discribe  the  progress  of  the  epidemic  when  n 
is  large.  There  is  no  apparent  reason  for  the  various  transition  rates  to 
have  the  same  proportionally  coefficient,  as  is  assumed  by  BLL(1979).  We 
modify  those  transition  rates,  given  by  Equation  (1.11),  and  assume  that  for 
te[0,<»),  and  r  *  1,  ....  m 


(6.2) 


RCX„(t),r)  =  -j 


-1 


V  Xn,rCt)(n-Xn{t)+Xn(0;))’  n  *  2> 


a  X  (t) 


n  *  «. 


Next,  we  derive  for  the  models  determined  by  Equation  (6.2)  the  desired 
approximations.  Let  V(b,y,t)  «  raaxU:  k  -  1 ,2, . . . ,£*=1 (b+q-l)-1UqsYt) , 
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b  =  1,  2,  . ..,  t,  Ye(0,“).  For  the  sake  of  completeness  we  present  the 
following  definition. 

Definition  6.1.  We  say  that  a  rva  W  is  negative  binomial  with  parameters 

a  =  1,  2,  ....  and  per  (0,1),  and  write  W'vNB(a.p)  if  P{W=k}=  a+k’1  Jpa  (l-p)K , 

l  a-i  j 

k  »  1,  2 . 

Note  that  by  BLL(1979)  p.  19B  we  obtain  that 

(6.3)  V(b,Y,t)^NB(b,e"n),  b  =  1,  2,  ....  t,  ye(0,-). 


he  show  next  that  for  all  te(0,°°),  and  r  =  1,  ....  m. 


m  -at 

(6.4)  fim  P  .  (t)  =  n  P(NB(X<b  (0),e  r  )=k 

n-w  »—  r=l  ’ 

-at 

(6.5)  £im  P  (t)  =  PiNBtX^  (0),e  r)=k), 

n»K»r  >r 


} ,  and 


Note  that  Conditions  (1.8)  and  (1.12),  with  6(k,x)  =  ark,  hold.  Consequently 
Statements  (6.4),  (6.5)  follow  by  Theorems  4.3,  5.3,  and  Statement  (6.3). 
Finally,  we  show  that  for  tc(0,°°),  r  =  1,  ...,  m 

V 

(6.6)  am  EX  =  X  (0)  (e  -1),  and  that 

n>r  °°>r 

2a  t  at 

(6.7)  am  Var{X  ft)}  =  l  ,.(0)  (e  r  -e  r  ) . 

n~°  n'r  *r 


In  particular  we  obtain  from  Statements  (6.6),  (6.7),  that  for  all  te(0,®) 

a  t 

(6.8)  am  EX  (t)  =  r  Xw  (0)  (e  r  -1),  and 

n-H»  ' 

2a  t  at 

(6.9)  am  Var(X  (t) }  =  I”  .X,  (0)  (e  r  -e  r  )  . 

n-*=  ’ 

To  prove  Statements  (6.6),  (6.7),  let  <T  *  Ir_iar*  Tben  for  all  tc(0,«)  and 
almost  all  nc(l,2, . . .}P(Xn(t)2k+Xn(0)}  s  Pfl^Cn"1  (n-q+1)  (Xn(0)+q-l) ]"lU  sot}  * 

s  P{V(Xw'0) ,a,t)2k},  k  ■  0,  1,  ...  .  Thus,  Conditions  (1.9),  (1.10),  hold  for 


all  t,  Be  (C ,") ,  with  Y(t)  =  V  (X^  (0)  ,a  ,t) .  Consequently  Statements  (6.6j, 
(6.7),  follow  from  Theorems  (4.6),  (5.3),  Statement  (6.3),  and  some  simple 
moments  evaluations  of  negative  binomial  tv's. 

Frequently  infectives  are  denied  the  freedom  of  movement ,  as  is  the  case 
when  they  are  hospitalised  immediately  after  they  contract  a  disease.  We 
describe  the  progress  of  this  epidemic  by  an  m-dimensional  stochastic  process 
determined  by  the  following  transition  rates: 


(6. 10)  RQCJt) 


Now,  we  derive  for  the  models  determined  by  Equation  (6.10)  the  desired  approx - 

.  a1 -i  a.-l  ,  m  . 

imations.  Let  S(a  ,a  )  =  / z  (1-s)  d/.,  and  (k.  , . . .  ,k  )  =  k !  [  n  k  l*1, 

r=  i  * 

kl'  •••»  k  :  for  the  sake  of  simplicity  we  obtain  these 

approximations  only  for  the  symmetric  case:  when  aj  *  a2  •• •  =  am  s  a •  Lemma 
(2.1)  it  follows  that  for  k  =  1,  2,  ... 


(6.1i)  P{1  =£  ,q«l,...,k)  = 


«  (X,  ..  (0)  +dr  (k )  ' 


«,q  q 


.  r=l  j 


dr(k) 


-  , 


,fx„(0).k-l 


re  k 


,  where  ,  ...,  ^  e{l,...,m},  and 


dr(k)  =  Iq3jI(^q=r)>  r  «  1,  ....  m. 

By  Lemma  2,2  we  conclude  that  the  interinfection  times:  T  T  are 

iid  exponential  rva's  with  a  me  equal  to  a”1  independent  of  the  infection 
causes:  ^ ,  ...  .  In  particular  note  that  X^t)  -  X^fO)  is  a  Poisson 


rva  with  parameter  ot ,  t« (0 ,«) .  Thus ,  for  t€ (0 ,«) ,  and  r  ■  1 ,  . . . , 


m 
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(6.12) 


(6.13) 


P~,k(t)  = 


-at,  . . s  (k)  ‘ 
e  (at)  —  H 

r=l 


X  (0)+k  -1 
0D,r'-  ‘  r 


'x.(o)+S(k)-iy 

(s  (k) ) : 

sQO  J. 

,  and 


-at , 


P«,k,r(t)  =  e  W  BC^^COJ.X.W-X.^CO))- 

•Iq=0(at)qS(X-.,r(0)+k*X»(0)-X«lr(0)^K^rl  = 

-1  rlT^r^-1  XcoCOJ-X^  (0)-l 

-  B  *(X  (0),X„(0)-X  (0)/A  z  *r  (l-z) 

*  1  ’  0  L 


,  „  ,k  -atz  -1 
(atz)  e  (k!) 


dz . 


Since,  Condition  (1.8)  holds  we  obtain  by  Theorem  4.3  and  Corollary  4.4  that 
pnk(t),  Pn,k,r ^  can  be  aPProximate  respectively  by  P^ft),  P  k  r(t),  given 
by  Equations  (6.12),  (6.13). 

Next,  note  that  for  te(0,°°)  and  all  net  1 ,2, . . .  }P{Xn  (t)>k+Xn  (0)  ^Sat), 

k  =  0,  1 .  Thus,  Conditions  (1.9),  (1.10)  hold  for  all  t,  te(0,»),  with 

Y(t)  being  a  Poisson  rva  with  mean  at,  te(0,“).  Consequently  by  Theorem  4.6 
and  Statement  (6.12)  we  obtain  for  te(0,“)  and  r  ■  1,  ...,  m,  that 

(6.14)  dim  EX  (t)  =>  atX^  (O)X"1  (0) ,  and  that 

n-H»  n,r  °° ,r 

(6. 15)  dim  Var{X  ft)}  =  atX  (0)X_1(0)* 

n-~  n,r  ,r 

♦  a2t2X->r(0)  (Xoo(0)-X<B>r(0))x;1  (0)  . 

Further,  by  Theorem  4.5  we  obtain  that  for  te(0,«) 

•dim  EX  (t)  «  at  and  dim  Var{X  (t)>  -  at. 

n+«  n-*» 


(6.16) 


Next,  we  consider  m-dimensional  simple  epidemic  models  given  by  the 
following  transition  rates 


(6.17)  RQ^CtJ.r)  = 


V’0Xn,r(t)(n'Xn(t)+Xn(O))0*  n  =  1,  2,  .... 

£ 

.a  X  (t)  ,  n  =  °°, 

(  r  «,r  ’  ’ 


where  ?,6<;(0,°°),  r  =  1,  ....  m. 

The  m-dimensional  simple  epidemic  models  determined  by  Equation  (6.17) 
generalize  univariate  simple  epidemic  models  used  and  motivated  by  Severo 
(1969).  Note  that  Conditions  (1.8)  and  (1.12),  with  A(k,r)  =  a^k,  hold. 
Consequently  by  Theorems  4.3,  5.3,  we  obtain  for  te(0,°°),  r  =  1,  . . . ,  m,  that 

(6.18)  £im  P  .  (t)  =  P{ 7k  ,  (X  (0)+q-l)'6U  <a  t), 

n»k»r  4q*lv  “>r  4  q  r  ’ 


and  that 
(6.19) 


Pn,klt)  ■  P„,k 

n-**>  r=ln-*<»  ’  r  ’ 

-6 


Let  W(t)  =  max{k:£  .  (X^OJ+q-l)  U  Sa  t).  Then  for  all  te(0,°°)  and  almost 
q— i.  q 

all  nc{l ,2 , . . . }P{Xn(t)sk+Xn(0) }  <  P{W(t)2k),  k  =  0,  1,  ...  .  Thus,  Conditions 
(1.9),  (1.10),  hold  for  all  t,  0e(O  ,<*>),  with  Y(t)  =  W(t) .  Consequently  the 
moments  approximations  hold  for  the  models  determined  by  Equation  (6.17). 

Finally  following  Gart  (1968),  (1972),  consider  a  population  of  sus- 
ceptibles  exposed  to  a  disease  and  partitioned  to  m  subpopulations.  Gart 
assumes  that  the  susceptibility  level  of  an  individual  varies  according  to 
his  membership  in  the  various  subpopulations.  With  small  adjustments  the 
process  X^Ct)  can  be  used  to  describe  the  progress  of  this  univariate  simple 
epidemic.  In  this  case  the  components  of  X^Ct)  represent  the  number  of  in- 
fectives  in  the  respective  subpopulations  at  time  t  measured  from  TQ.  Further, 
all  the  results  obtained  in  this  paper  carry  through  for  these  univariate  models. 
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